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Abstract: We compare the relic density of neutralino dark matter within the minimal 
supergravity model (mSUGRA) using four different public codes for super symmetric spec- 
tra evaluation. While the predictions for the relic density of neutralinos are rather stable 
in most of the mSUGRA space, it is in the most physically interesting regions that large 
discrepancies can be observed, in particular the focus point, large tan/3 and coannihilation 
regions. 
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1. INTRODUCTION 

One of the most stringent constraints on supersymmetric models with R-parity conservation 
arises from the upper hmit on the rehc density of dark matter. This is particularly true 
with the recent precise measurements of the cosmological parameters realised by WMAP. 
It is therefore crucial to quantify the theoretical uncertainties that enter the calculation 
of the relic density of the lightest supersymmetric particle (LSP) and to see how they 
reflect on the allowed parameter space. We do not attempt to answer this question fully 
here. We will only consider one aspect: the uncertainty introduced by the calculation 
of the weak scale SUSY parameters using renormalization group equations (RGE) within 
the context of the mSUGRA model. As a measure of the theoretical uncertainty on the 
mSUGRA parameters, we use the four public state-of-the-art RGE codes: Isajet7.69 Q, 
SOFTSUSYl.8.3 [|], SPHEN02.20 [| and Suspect2.2 [§, link them to micrOMEGAsl . 2 §] 
and compare estimates for the relic density. At this point no attempt is made to estimate 
the uncertainties that could arise directly in the calculation of the relic density itself. 

2. RGE CODES AND RELIC DENSITY CALCULATIONS 

A detailed study of theoretical uncertainties on the supersymmetric spectra as obtained 
by RGE codes was presented in |^. It was shown that differences in masses less than a 
few percent are usually found, although some corners of parameter space are still difficult 
to tackle and can display much larger differences. The discrepancies can be traced back 
to the level of approximation used in the weak-scale boundary conditions. The large tan (3 
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region and the focus point region (large Mq) are still subject to large theoretical errors. 
Both of these regions are precisely where one can find cosmologically interesting values for 
the relic density, O/i^ < .128. In the focus point region, the LSP is mainly a Higgsino 
and annihilates efficiently into gauge bosons. At large tan f3, even rather heavy neutralinos 
can annihilate into bb pairs via s-channel exchange of a heavy Higgs. The coannihilation 
region where the Next-to-Lightest super symmetric particle (NLSP) is nearly degenerate in 
mass with the LSP, is another cosmologically relevant region. Although it is a priori not 
difficult to handle by the RGE codes, the value of the relic density depends sensitively on 
the mass difference between the NLSP and the LSP and even shifts of GeV can cause 
large shifts in the relic density. The other cosmologically viable mSUGRA region, the bulk 
region, shows a much smaller induced sensitivity upon the MSSM mass spectrum. 

The link between micrOMEGAsl .2 and the RGE codes is done within the spirit of the 
SUSY Les Houches Accord : common input values are chosen and pole masses, mixing 
matrices, the parameter and the trilinear couplings are calculated by the RGE codes. All 
parameters are read by micrOMEGAsl . 2. The annihilation cross-sections are then evaluated 
at tree-level. Important radiative corrections to the Higgs widths and in particular the Amb 
correction are taken into account. 

3. RESULTS 

For the numerical results as default values we have fixed rrit = 175 GeV, as(M^)*^'^ = .1172 
and mh{mb)^'^^ = 4.16 GeV. This corresponds to mh{Mz)^^ = 2.83 GeV. We concentrate 
on the three regions where the relic density is within the WMAP range and where poten- 
tially large discrepancies can be observed: the focus point region, the large tan /? region 
and the coannihilation region. 

3.1 Coannihilation 

Mo = 150 GeV, = 0, tan /? = 10, > 

The small M1/2 region corresponds to the so-called bulk region where the bino-LSP anni- 
hilates into lepton pairs via s-channel Z or Higgs exchange or t-channel slepton exchange. 
Here one finds very good agreement between the values of O/i'^using the different RGE 
codes (see Fig. ||a) since the predicted values for slepton and neutralino masses are in good 
agreement (within a few GeV). The exact position of the Z pole (corresponding to the big 
dip in Qh'^) is slightly shifted for SPHEN02.20 but the range of values of Mi^2 which 
nh^ < .128 are basically identical. Note that the Z pole region is ruled out by the LEP 
constraints on neutralinos within the context of mSUGRA models. 

As one moves up in M1/2, one reaches the so-called coannihilation region where the f 
is the NLSP and is nearly degenerate with the neutralino, as in Fig. |T|b. Coannihilation 
with the f, and to a lesser extent the selectron and smuon, brings the relic density in the 
desired range. For a given value of M1/2, differences between the codes can reach a factor 2, 
the largest differences are found between SPHEN02.20 and SOFTSUSYl . 8 . 3. However very 
good agreement is found between all codes when the relic density is plotted as a function of 
the mass difference between the LSP and the NLSP (here the f). All codes obtain values 
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Figure 1: a) nh'^ for Mq = 150 GeV, tan^ = 10, = 0, ^ > for SOFTSUSYl . 8 . 3 (full), 
Isajet7.69 (dashed), Suspect2.2 (dash-dotted), and SPHEN02.20 (dotted). At large M1/2, 
Isajet7.69 and Suspect2.2 give nearly identical results, b) flh^ vs mf^ — m-^o for the same 
set of parameters as a). 



of Sl/i^compatible with WMAP for mass differences — m^o ~ 4 GeV (at the extreme 
left of Fig.|l|b), even though the corresponding value of the neutralino mass can differ. The 
value of Ml /2 for which the relic density becomes compatible with WMAP varies from 670 
GeV (SPHEN02.20) to 790 GeV (SDFTSUSYl . 8 . 3), a 12% difference on M1/2. 

3.2 Focus point 

Ml/2 = 300 GeV, Aq = 0, tan/3 = 10, /i > 

In addition to the small Mq (bulk/coannihilation region) where annihilation into leptons 
is important, the cosmologically relevant region is found at values of Mq well above ITeV. 
As one approaches the region where electroweak symmetry breaking is forbidden, the 
parameter approaches zero. This means that the LSP is mainly Higgsino. This LSP can 
then annihilate very efficiently into gauge bosons (WW/ZZ) and to a lesser extent into Zh. 
The parameter fi is however very sensitive to the top Yukawa coupling, ht (which is also 
reflected in a sensitivity to the value of the top quark mass) and huge differences between 
codes were observed Q. The impact on the relic density and on the exclusion region is 
likewise very significant. 

As can be seen in Fig. ^, all codes agree very well for Mq < ITeV but as one gets to 
large values of Mq, more than one order of magnitude differences in can be found. 

For rrit = 175 GeV, only Isajet finds a large drop in the /i parameter as one moves to 
Mq ~ 3000 GeV, this is when drops below the upper limit from WMAP. The other 
codes do not find this drop in /i and do not obtain a cosmologically interesting region for 
Mq < 4000 GeV. These large differences between codes however are just a reflection of the 
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Figure 2: a) fi/i^ as a function of Mq for A/1/2 = 300 GeV, tan (3 ^ 10, Aq ~ and /i > and mt — 
175. Same labels as in Fig. 1. b) Dependence of the relic density on rrit for SOFTSUSYl . 8 . 3 (full) 
and SPHEN02 . 20 (dash) . 



sensitivity to the top Yukawa, ht{MsusY) which is proportional to rrit. We show in Fig. 
2b, the variation of Q.h'^ with rrit using SOFTSUSYl . 8 . 3 and SPHEN02.20 for Mq = 3000 
GeV. The value Vth^ = .128 found in Isajet7.69 for nit = 175 GeV can be reproduced in 
SOFTSUSYl. 8. 3 (SPHENO) by changing the input to nit = 172.2(172.5) GeV. 



3.3 Large tan/3 

mi/2 = 1500 GeV, Aq = 0,tan^ = 52 > 

At large tan (3 the new feature is the annihilation of neutralinos into bb via heavy Higgs 
exchange. With the current version of the RGE codes, this is observed only for very large 
values of tan (3. The crucial parameter here is M_4/2m^o which must be close to unity to 
provide sufficient annihilation of neutralinos. Large differences in the value of Ma between 
the different RGE codes occur because of the sensitivity of the RGE to the bottom Yukawa 
as well as from taking into account higher loop effects. 

As Fig. 3a shows, all 4 programs predict a large drop in the relic density when the 
neutralino mass gets close to Ma/'^ although this drop occurs at much lower values of 
Ml/2 for SPHENO, M1/2 1250 GeV than for Isajet7.69 , M1/2 ~ 1750 GeV. However, 
here again the results are very sensitive to the input parameters, in this case the value of 
the b-quark mass. For M1/2 = 1300 GeV, we find an order of magnitude shift in Q.h? for 
nih{nih) =4-4.4 GeV with the program SOFTSUSYl . 8 . 3. By a slight shift of the b-quark 
mass we can find perfect agreement between SPHEN02.20 and SOFTSUSYl . 8 . 3, as shown 
in Fig. 3b. 
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Figure 3: a) flh^ as a function of A/1/2 for nit ~ 175, tan/3 = 52, Aq — and /i > 0. Same 
labels as in Fig. 1. b) Dependence of the relic density on mi,{mi,) for SOFTSUSYl . 8 . 3 (full) and 
SPHEN02.20 (dash). 



4. CONCLUSION 

While the predictions for the relic density of neutralinos are rather stable in most of the 
mSUGRA space, it is in the most physically interesting regions that large discrepancies 
can be observed, in particular the focus point, large tan/3 and coannihilation regions. It is 
however reassuring to find that with the newer versions of the codes, the discrepancies in 
the sparticle spectra tend to be reduced. More details on the theoretical uncertainties in 
the evaluation of the relic density arising from the standard model parameters, as,mb, mt, 
used as input in a RGE code can be found in Q . 
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